INVERSE SCATTERING TRANSFORM 

FOR THE DEGASPERIS PROCESI EQUATION: 

A RIEMANN HILBERT APPROACH 
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Abstract. We present the inverse scattering transform approach to the Cauchy problem 
on the line for the Degasperis-Procesi equation 
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in the form of an associated Riemann-Hilbert problem. This approach allows us to give 
a representation of the solution to the Cauchy problem, which can be efficiently used in 
studying its long-time behavior. 



xjy . 1. Introduction 

In this paper we present the inverse scattering approach, based on an appropriate Riemann- 
Hilbert problem formulation, for the initial value problem for the Degasperis-Procesi (DP) 
equation [TTJ, [16] 

u t - u txx + 3uju x + Auu x = 3u x u xx + uu xxx , -co < x < oo, t > 0, (1.1) 

u(x,0) = u {x), (1.2) 



> 

in 

o\ 

Q\ . where wisa positive parameter. The DP equation arises as a model equation describing the 

shallow-water approximation in inviscid hydrodynamics in the so-called "moderate amplitude 
regime" : introducing two small parameters, the wave-amplitude parameter e (characterizing 
the smallness of the wave amplitude) and the long-wave parameter 5 (characterizing the 
smallness of the typical wavelength with respect to the water depth), in this regime we 
assume that 5 <C 1 and e ~ 5. This regime can be characterized as to be more nonlinear 
than dispersive, which, in particular, allows "wave breaking". This is in contrast with the 
so-called "shallow water regime" (<5<1 and e ~ <5 2 ), where nonlinearity and dispersion are so 
balanced that the solution of the initial value problem for the associated nonlinear equation 
(the Korteweg-de Vries equation) exists globally for all times, for all "nice" (sufficiently 
decaying and smooth) initial data. 

Among the models of moderate amplitude regime, only two are integrable (admitting a 
bi-Hamiltonian structure and a Lax pair representation): they are the Camassa-Holm (CH) 
equation and the DP equation. Also, they are the only two integrable equations from the 
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2 A. BOUTET DE MONVEL AND D. SHEPELSKY 

"6-family" of equations 

ut ~ utxx + buu x + (b + l)uu x = bu x u xx + uu xxx . 

The CH and DP equations correspond to b = 2 and 6 = 3, respectively. 

The analysis of the CH equation by using the inverse scattering approach were started in 
|12 |, I15 |, [20] . A version of the inverse scattering method based on the Riemann-Hilbert (RH) 
factorization problem was proposed in [5j [7J (another Riemann-Hilbert formulation of the 
inverse scattering transform is presented in [13]). The RH approach has proved its efficiency 
in the study of the long-time behavior of solutions of both initial value problems [USE] and 
initial boundary value problems [8j. 

In the present paper we develop the Riemann-Hilbert (RH) approach to the DP equation, 
following the main ideas developed in [7J. 

A major difference between the implementations of the Riemann-Hilbert method to the 
CH equation and the DP equation is that in the latter case, the spatial equation of the 
Lax pair is of the third order, which implies that in the matrix form one has to deal with 
3x3 matrix-valued equations, while in the case of the CH equation, they have a 2 x 2 
matrix structure, as in the cases of the most known integrable equations (KdV, modified 
KdV, nonlinear Schrodinger, sine-Gordon, etc.) Hence, the construction and analysis of the 
associated RH problem become considerably more complicated. 

In a recent paper [II], the inverse scattering method for the DP equation based on a 3 x 3 
matrix RH problem in the spectral /c-plane is proposed, where the solution of the DP equation 
is extracted from the large-A: behavior of the solution of the RH problem. Also, the dressing 
procedure is given for constructing iV-soliton solutions, which is illustrated, particularly, by 
the explicit construction (in parametric form) of the 1-soliton solutions. 

In our approach, we propose a different Riemann-Hilbert problem and give a different 
representation of the solution u(x, t) of the initial value problem (|1.1|) in terms of a solution 
of this problem evaluated at a distinguished finite point of the plane of the spectral parameter. 
Remarkably, the formulae for u(x,t) have the same structure as the parametric formulae for 
pure multisoliton solutions obtained in [21] . 



2. Lax pair and eigenfunctions 

We assume that the initial function uq(x) in f)1.2|) is sufficiently smooth and decay fast as 
\x\ — > oo. Letting u(x,t) be the solution of the DP equation, we introduce the "momentum" 
variable 

m = m(x, t) = u — u xx . 
It is known that, similarly to the case of the CH equation (see, e.g., [32]), the condition 
m(x, 0) + u > for all x provides the existence of m(x, t) for all t; moreover, m(x, t) + uj > 
for all x G R and all t > 0. This justifies the form 

{^T+^) t = - (u^+ZP) x (2.1) 

of the DP equation, which will be used in our constructions below. 

The linear dispersion parameter uj can be scaled out to u> = 1. Hence, in what follows, for 
simplicity, we assume u> = 1. 
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2.1. Lax pairs. 

Scalar Lax pair. The DP equation admits a Lax representation: this equation is actually the 
compatibility condition of two linear equations [16] 

i>x ~ tpxxx = -z 3 (m + 1)^, (2.2a) 

ip t =(u x g-J-0 -utf} x + — il} xx , (2.2b) 

where z is a spectral parameter, ip = ip{x, t, z), and c is an arbitrary constant. In what follows 
we will see that it is convenient to choose c = §. 

1st matrix form. In order to control the behavior of solutions to (|2.2[) as functions of the 
spectral parameter z (which is crucial for the Riemann-Hilbert method), it is convenient to 
rewrite the Lax pair in matrix form. Introducing $ = $(x, t, z) by 

i> 

transforms the scalar Lax pair (|2.2p into a Lax pair 

$ x = U3>, (2.3a) 

$t = F$, (2.3b) 

where 

U(x,t,z) = I Oil, (2.3c) 
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V(x,t,z)=[ u + l -Bj£ -u I (2.3d) 

with 

q = q(x,t) = (m+l) 1/3 . (2.3e) 

From now on, (|2.3|) will be seen as a 3 x 3 matrix-valued linear system: a "matrix" solution 
of (|2.3|) is a collection $ = (<1>W $( 2 ^ <&( 3 )) of three "vector" solutions $( J ). Now we notice 
that if c = 3, then V is traceless. Thus in this case the determinant of a matrix solution to 
the equation $t = V<& is independent of t. The analogous property of <& x = [7<3? is obvious. 

The coefficient matrices U and V in (I2.3P have singularities (in the extended complex z- 
plane) at z = and at z = oo. In order to control the large-z behavior of solutions of (J2.3J) . 
we follow a strategy similar to that adopted for the CH equation [5j [7] . We transform (|2.3p 
in such a way that: 

i) the leading terms for z — > oo in the Lax equations be diagonal, whereas the terms of 

order 0(1) be off-diagonal; 
ii) all lower order terms (including those of order 0(1)) vanish as |rc| — > oo. 

It is instructive to perform this transformation in two steps: 
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(i) First, transform (|2.3p into a system where the leading terms are represented as products 

of (x, £)-independent and (x, independent factors, 
(ii) Second, diagonalize the (x, £)-independent factors. 

2nd matrix form. For the first step, introducing $ = <t(x,t,z) by 

where 

D(x, t) = diag{g _1 (x, t), 1, q(x, t)} 
transforms (|2.3p into a new Lax pair 

§ x = U$, (2.4a) 

<l t = T/i, (2.4b) 

where 

/0 1 o\ /f o 

U{x,t,z) =q{x,t) 1+0 

V 3 i o/ \o |- 

= g(x,t)C/ 00 (z) + L>( 1 )(x,t) (2.4c) 

and 

/0 1 N 
F(x,i,2) = -uq 1 | + 




= -u(x, t)q(x, QUoo (z) + V^ (z) + V^i) (x, i) + ^T/( 2 ) (x, t) . (2.4d) 

Main matrix form. For the second step, it is important that the commutator of C/qo and Yx 
vanishes identically, i.e., [t/"oo)Kx>] — 0) which allows simultaneous diagonalization of U^ and 
Voo- Indeed, we have 

p- 1 (z)U OD (z)P(z)=A(z), (2.5a) 

P- 1 (z)V 00 (z)P(z) = ±I + A-\ z ) = A(z), (2.5b) 



z 3 ' 



where / is the identity 3x3 matrix, 



%(z) 

A(z) = ( X 2 (z) | , (2.6a) 

X 3 (z) / 

1 1 1 
P(z) = | Ai(z) \ 2 {z) X 3 (z) | , (2.6b) 

\l(z) \%(z) A§(*), 
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and 

/ {SX^zyi)- 1 \ (\ftz)-l Ai(*)l\ 

P (z) = (SXziz)-!)- 1 Al(«)-1A 2 (^)1 . (2.6c) 

V (SAsM-l)- 1 / \A|( Z )-1 X 3 (z) l) 

Here Aj(z), jf = 1, 2, 3 are the solutions of the algebraic equation 

A 3 - A = z 3 , (2.7) 

so that Xj(z) ~ CjZ as z — > oo, where the Cj's are the cube roots of unity. 
Thus, introducing $ = &(x,t,z) by 

6 = P- 1 ® 

leads to another Lax pair 

$ x - qA(z)$ = U$, (2.8a) 

$ t + (uqA(z) - A(z))$ = V$, (2.8b) 

where 

U{x,t,z) =p- 1 (z)U {1) {x,t)P(z), (2.8c) 

V(x,t,z) = p-\z) (v m (x,t) + -iy( 2 )(x,t)j P(z). 



(2.8d) 



Commutator form. Notice that U(x,t,z) = 0(1) and V(x,t,z) = 0(1) as z -^ oo since 
LA 1 ) and y' 1 ) are lower triangular matrices. Moreover, it can be checked directly that the 
diagonal entries of U(x,t,z) and V(x,t,z) are of order 0(1/ z). This latter fact is important 
to establish the large- z behavior of $ pQ. 

On the other hand, U(x,t,z) = o(l) and V(x,t,z) = o(l) as |x| — > oo, which suggests 
introducing a 3 x 3 diagonal function Q(x, t, z) that solves the system 

Qx = qA(z), (2.9a) 

Qt = -uqA(z) + A(z), (2.9b) 

as follows: 

Q(x,t,z)=y(x,t)A(z)+tA(z) (2.10) 

with 

y(x,t) = x- / (?(£,*) -l)d£. (2.11) 

J X 

The fact that the equations in (|2.9p are consistent follows directly from the DP equation in 
the form f)2. lj) : % = —(uq) x . The normalization of Q(x,t,z) is chosen in such a way that 

Q(x,t,z) ~xA(z)+tA(z) as i ^ +00. (2.12) 

The role of Q(x,t,z) in the construction of the Riemanm-Hilbert problem is to catch the 
large- z behavior of dedicated solutions of the Lax pair equations (|2.8p . Indeed, introducing 
the 3x3 matrix-valued function M = M(x, t, z) by 

M = l>e~ 9 
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reduces (12. 8ft to the system 

M X -[Q X ,M] = UM, 
M t -[Q t ,M]=VM. 

2.2. Eigenfunctions. 



(2.13a) 
(2.13b) 



Fredholm integral equations. Particular solutions of (|2.13p having well-controlled properties 
as functions of the spectral parameter z can be constructed as solutions of the Fredholm 
integral equation (cf. [1]) 

M(x,t,z) = 

<x,t) 



1 + 



^Q(x,t,z)-Q(£,T,z) 



(x*,t*) 



UM((, t, z)d£ + VM((, t, z)dr e 



-Q(x,t,z)+Q(£,T,z) 



(2.14) 



where the initial points of integration (x*,t*) can be chosen differently for different matrix 
entries of the equation: Q being diagonal, (J2.14H must be seen as the collection of scalar 
integral equations (1 < j, I < 3) 



Mji(x,t,z) 

r(x,t) 



^^''• 2 )- Q ^(«' r ' 2 )((L>M) ji (^r,z)de + (yM) j i(C,r,z)dr)e- Q " (:c ' i ' 2)+Q " ( «' T ' 2) . 
'(*Jj.*Ji) 

Choosing the (x^^t*^) appropriately allows to obtain eigenfunctions that can be used in the 
construction of the Riemann-Hilbert problems associated with the initial value problems [7] 
as well as the initial boundary value problems [8] . 



->— T 



t ----?- 







x £ x £ 

Figure 1. Paths of integration. Left: (x*,t*) = (— oo,i). Right: {x*,t*) = (+oo,t). 

In particular, for the Cauchy problem considered in the present paper, it is reasonable to 
choose these points to be (— oo, 0) or (+oo, 0) thus reducing the integration in (|2.14|) to paths 
parallel to the x-axis (see Figure [Q): 



^Q(x,t,z)-Q(£,t,z)i 



-Q(x,t,z)+Q(£,t,z) 



M(x, t,z)=I+ I e V\*M-WM [UM(£, t, ^WJ+viwj^ 

(±)oo 



or, in view of (|2.10p . 
M(x,t,z) = I + 



(±)oo 



e- /J(<KC,*)-i)dCA« p M{ ^ tj 3)] e /jfo(C,t)-i)dCA(z) dC> 



(2.15) 



(2.16) 
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Since q — 1 > 0, the domains (in the complex z-plane) where the exponential factors in (12.15P 
are bounded are determined by the sign of Re\j(z) — ReXi(z). 

A new spectral parameter. It is convenient to introduce a new spectral parameter k (see [14] 
and also [TT]) such that 

1 / 1\ 1/3 
z(k) = -=k 1 + T7? > (2-17) 



(2.18) 



and fixed by the condition z(k) ~ -4= k as fc — > oo. We fix Xj = Xj(k) as follows 

Xn(k) = —= \ uj^k -\ r- | where ui = e~3~. 




Figure 2. Rays /„, domains f2„ and points x w in the fc-plane 
The set E = {A; | Re Aj(fe) = Re A; (A;) for some j ^ 1} consists of six rays 



i(u-i) 



v = 1,... ,6 



l v = M+e^ 
dividing the A>plane into six sectors 

{TT TT 1 

k ^(^- 1) < arg A; < -i/|, z/ = l,...,6. 

In order to have a (matrix- valued) solution to (J2.16P to be analytic in C\S, the initial points 
of integration oo,-/ are specified for each matrix entry (j, I), 1 < j, I < 3 as follows: 

+oo, if KeXj(z) > ReXi(z), 
— oo, if ReXj(z) < ReA;(z). 

That means that we consider the system of Fredholm integral equations 

M jl (x,t,z)=I jl + f X e-^'^i'i^^-^ dC [(t>M) j7 ^,t^)]e Ai ^i'i^( C '*)- 1 ) dC de. (2.20) 



oo n 



(2.19) 
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Proposition 2.1 (analyticity). Let M(x,t,k) be the solution of the system of equations 
(I2.20p . where the limits of integration are chosen according to (|2.19p . Then, 

(i) M is piecewise meromorphic with respect to £, as function of the spectral parameter k. 
(ii) M(x, t, k) -^ I as k ->■ oo. 

Proof. The proof follows the same lines as in pQ. Notice that in order to have (ii) it is 
important that the diagonal part of U vanish as k — > oo. □ 

Proposition 2.2 (symmetries). M(x,t,k) satisfies the symmetry relations: 

(51) F 1 M(x,t,k)T 1 = M(x,t,k) where F t = (iooY 

(52) T 2 M(x,t,kuj 2 )T2 = M(x, t, k) where F 2 = ( o l Y 

(53) T 3 M(x,t,kuj)r 3 = M(x, t, k) where r 3 = ( o o l Y 

Hence, the values of M at fc and at wA; are related by 

M{x,t,ku) = C- 1 M(x,t,k))C where C=(ioo). 

If Xi(k) = Xj(k), i 7^ j for certain values of the spectral parameter k, then P at these values 
becomes degenerate (see (|2.6cp ). which in turn leads to singularities for U and, consequently, 
for $ and M. These values are x„ = e~3~V J ~ ', v = 1,...,6. Taking into account the 
symmetries described in Proposition 12,21 leads to the following 

Proposition 2.3 (singularities). The limiting values of M(x,t,k) as k approaches one of 
the points x v = e"'^ -1 ^, u = 1, . . . ,6 have pole singularities with leading terms of a specific 
matrix structure. 

(i) At k = >c\ = 1, as k = 1 ± ie — >■ 1, e > 0, 

M~ | -i -1 -J 1 | (I n+ ! • 12.21) 

where a_ = — cEjT and /3_ = — f3 + . 

/ \ — — — 

(ii) At k = X2 = e 3 , as k = e 3 ± ie — !• e 3 , £ > 0, 

I \ (J3 ± 
M ~ ^^ (1 1 1 d± | , (2.22) 

1 -1 -1/ \ a ±/ 

where a- = —cja+ and /3- = —loj3 + . 

M has similar leading terms at the other polar singularities k%, . . . , kq in accordance with the 
symmetry conditions stated in Proposition \2.S[ 
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Indeed, consider, for example, the behavior of M at k = 1. We have Ai(l) = A2(l) 
— 1/\/3 and \s(l) = 2/\/3- Moreover, as k — >• 1, we have 

3A? - 1 = -2\/3(& - 1) + O ((Jb - l) 2 ) , 
3A 2 . - 1 = 2V3(/fc - 1) + O ((jfe - l) 2 ) , 
3A§ - 1 = 3. 

Consequently, at k = 1, 

P_1 ( fc ) = r^r ! - J - J -' II ° /'-' n I +°( J -) 




with some pj . 

Now recall that $ = P _1 $, where, on one hand, 6 satisfies a differential equation whose 
coefficients are regular at k = 1, and on the other hand, $ satisfies the boundary condition 
(see (|2.16j) with signs chosen according to (|2,19p ) 

4> ~ P ( k ) e y^t)K{k)+tA(k) ; x ^ +00 

for all fc such that Xj(k) ^ Xi(k) for all j ^ I. It follows that when k approaches 1 from the 
either side of 1%, the first two columns of $ coincide. Consequently, for the limiting values 
of $ (and thus M) we obtain (I2.2ip while the property a_ = — cEjT and /?_ = — j3+ follows 
from the symmetry (SI) of Proposition 12.21 Similarly for the other points k v = e"3"(" ', 
i/ = 2, . . . , 6. 

Remark 2.4. While the set of singularities of M(x, t, k) in the open domain 

fi = C\E = niU'"Ufi 6 

can be empty (for instance, this happens for all sufficiently small "potentials" u(x,t)), the 
singularities described in Proposition 12.31 are generic. This should be compared with |14j . 
where the solutions of a system of integral equations similar to (12. 16ft are combined into a 
Riemann-Hilbert problem under the assumption that they have no singularities. 

3. RlEMANN-HlLBERT PROBLEM 
3.1. Jump conditions. 

1st RH problem. For k on the boundary between the adjacent domains Q, v , the limiting values 
of M, being the solutions of a system of differential equations f|2. 13[) must be related by a 
matrix independent of x and t. Supplying the rays l v with the orientation, see Figure [21 we 
can write for the limiting values of M: 

M+(x,t,k) = M^(x,t,k)e Q( - x ' t ^S(k)e- Q ^ xAk \ kehU---Ul 6 . (3.1) 

Considering (I3.ip at t = we see that the conjugating matrix S(k) is, in fact, determined by 
u(x, 0), i.e., by the initial data for the Cauchy problem fll.lj) — <11.2j) . via the solutions M(x, 0, k) 
of the system (I2.16P whose coefficients are determined by u{x, 0). Thus the relation (13. If) can 
be considered as a "pre-Riemann-Hilbert problem" associated with (jl.ip : the data are S(k), 
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and we seek for a piecewise meromorphic function M satisfying (|3.1i for all x and t, in the 
hope that one can further extract the solution u(x,t) to (jl.ip from M(x,t,k). 

The conjugating matrix S(k) has a particular matrix structure, cf. [14]. Indeed, the 
integral equations (12. 16H allows studying the limiting values of M as x — > ±oo. Set t = 
and consider, for example, the limiting values of M±(x, 0, k) for k £ l\ = M + . The triangular 
structure of integration in 112. 16j) implies that 



Q(xM)M + (x,0,k)e Q( - x '°'V >S?(k)=\Q 1 (3.2) 



x— >-+oo 




and 



-Q(x,o,fc) M _( X)0;fe ) e Q(^o,fc) > £-(£) = I r_(Ar) 1 ] . (3.3) 




Then the symmetry (SI) from Proposition 12.21 implies that r~(k) = r+(k), and the jump 
matrix S(k) from (|3.ip for k G l\ takes the form 

S(k) = (S^{k)Y l S+{k)= \7(k) 1 | ( 1 0|, (3.4) 

where r(/e) := — r+(k). Similarly, for k G Z4 





-r(fe) 0\ 

5(fc) = (5r(fc)) -1 ^(fe) = [ r(fc) 1 I [ 1 , (3.5) 

whereas the construction of S(k) for k £ l„ with v ^ 1,4 follows from the symmetries of 
Proposition 12.21 Thus, similarly to the cases of, say, the KdV equation or the Camassa-Holm 
equation, the jump matrix on the whole contour is determined by a scalar function — the 
reflection coefficient r{k) given for l;£l. In turn, as it follows from (|3.2j) and (J2.16P for 
t = 0, the reflection coefficient is determined by the initial condition uo(x). 

2nd RH problem. The dependence of the matrix e^Se - ^ relating M + and M_ in (|3.ip on 
the parameters x and t suggests introducing the parameter y = y(x,t) given by (|2.1ip and 
thus rewriting (|3.ip as 

M+{y,t,k)=M-(y, t, k)e yA( -^ +tA ^ S{k)e- yMk) - tA( - k) (3.6) 

so that 

M(x,t,k) = M(y(x,t),t,k), 
which provides explicit dependence of the relating matrix on the parameters (y,t). 

3rd RH problem. On the other hand, the particular matrix structure of the singularities at 
k = x v = e"3"^ _ ' suggests introducing a row vector- valued 1x3 Riemann-Hilbert problem 

H+(y,t,k) = V-(y,t,k)e-y A ^- tA ^S(k) e y A ^+ tA ^ (3.7a) 
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having the same jump conditions as in the matrix case but normalized at k = oo by the 
condition 

fi(y,t,k) = (l 1 l)+o(l) as/c^oo. (3.7b) 

The transition from the 3x3 matrix RH problem to the 1x3 row vector RH problem can be 
viewed as the multiplication of the former by the constant row vector (111) from the left, 
which suppresses the singularities at k = >c v = e"^ -1 ). 

3.2. Residue conditions. In order to complete the formulation of the RH problem, one 
has to complete the jump condition (|3.7ap and the normalization condition (|3.7bp with the 
residue conditions at the possible poles of fi(y, t, k) or M(y, t,k) in C \ E where £ = U® =1 l u . 
Generically (see [T]), there are at most a finite number of poles lying in C \ £, each of 
them being simple, with residue conditions of a special matrix form: distinct columns of M 
(distinct entries of //) have distinct poles, and if k n is a pole, then 

Res fc=fcn n(y, t, k) = M (y, t, k n )^ k ^ +tA ^ v n ^y K ^- tA ^ , (3.7c) 

where the 3x3 matrix v n has only one non-zero entry at a position depending on the sector 
of C \ £ containing k n . For example, if k n £ fii, a non-zero entry of v n can be either (v n )i2 
or {v n )23- Then the positions (as well as the values) of the non-zero entries of poles in other 
sectors of C \ £ are determined by the symmetries (S1)-(S3). 

Similarly to the jump matrix S(k), the residue conditions are determined by the initial 
values uo(x) (putting t = in fl3.7c|) ). 

A closer look at the (y, independence of the exponential in the residue condition (|3.7c|) 
reveals the following. Suppose that the non-zero entry of v n is (v n )ji, and let 

v = Xj - A/. 

Then this exponential has the form (see f)2.5|) ) 

giXt-Xti+HAj-Ai) = e ^- Xl \ y -m) = e v (v-T^) } (3. 8 ) 



The last equality in (|3.8p follows from the fact (see the definition (12.7P of Xj) that 3AjA; = 
1 - (Xj - Xi) 2 for all j + I. 

An analogy with the Camassa-Holm equation suggests making the conjecture that the 
actual positions of the poles associated with a global solution of the initial-value problem 
(jl.ip with smooth, decaying initial data uo(x) are those for which the associated value of v 
is real and, moreover, < v < 1. 

Due to the symmetry relations, it is enough to consider the location of the poles in one 
sector of C \ E, say, in U\. 

The exponential factor for (1^)12 contains 

Ai(fc) - X 2 (k) = —= ( uk-\ --uj 2 k 



^/3 \ ook u 2 k 

OJ-LO 2 ( , 1\ . /, 1 



■ji V kj-v *.'■ (3 ' 9) 
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It follows that the l.h.s. of (13.9H is real either for {k \ Hek = 0} or for {k \ \k\ = 1}. The 
real axis does not intersect with Qi whereas for k = e lip with < ip < vr/6, the associated v 
satisfies the inequalities < \v\ < 1. 

The exponential factor for (v n )23 contains 

*.(*) - *.(*) = ^ (^ + ;?*-*- i) " '("^ =s) ■ ( 3 - 10 ) 

Thus the admissible arc for the location of the poles in f^i for this entry (23) is {k \ k = 
e lif , 7r/6 < (p < 7r/3}. Notice that the symmetry relations provide that the associated residue 
conditions in O3 are those for the (12) entry. Similarly for the other sectors £l u . 

The poles with the residue conditions described above are associated with the soliton 
long-time behavior of the solution of the Cauchy problem (|l.ip - (|1.2p . the velocity of all 
solitons being greater than 3 (since |t>| < 1 -- the velocity c of a soliton being related to 
v as c = 3/(1 — v 2 ), see (|3.8p ). On the other hand, the residue conditions with poles at 
{k I Re/c = 0} give rise to the "loop solitons" [22J : in this case, 

v = l(ip--J =-(p+-j, p£K 

and thus \v\ > 2; so they move, as opposite to the smooth solitons, in the negative direc- 
tion. These solutions are not classical ones: in the (y, t) scale, they are given, similarly to 
the solitons, in an univalent way, but the transition to the original (x, t) scale makes them 
multivalued. 

3.3. Solution of the Cauchy problem in terms of the solution of the RHP. 

Assumption. In what follows we assume that the RH problem consisting in finding a piece- 
wise, vector- valued function [i satisfying the jump condition (I3.7ah . the normalization condi- 
tion (|3.7bp . and the appropriate residue conditions has a unique solution. 

In order to obtain u(x,t) from its solution /j,(y,t,k), it turns out that it is convenient to 
evaluate /i(y,t,k) at distinguished finite points in the fc-plane, more precisely, at the points 

Ku = e T + e - {u ~ 1 >, I/ = l,...,6 
characterized by the property z(k, v ) = 0. 

Last Lax pair. Coming back to the system (|2.3p . let us introduce <$?(°' = &°'(x, t, z) by 

This reduces (|2.3p to 

¥V-A(z)^=U^&°\ (3.11a) 

^ 0) -A(z)^=V^¥°\ (3.11b) 
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where 

3\'i(z)-l ° 

m°\x,t,z)=z 3 m(x,t)\ ^^ | [ 1 1 1 ] , (3.11c 

0~3 




u -u ] P(z). (3.11d) 

■u^ — z 3 u(m + 1) 

Now notice that U^°'(x,t, z)\ = 0. Therefore, introducing M^> = M^°'(x,t, z) by 

M (o) = $(o) e -*A-tA 

and determining M^> as the solution of a system of integral equations similar to the system 
(|2.20|) determining M, we have: 

M^(x,t,z)\ z=Q ^I. 

On the other hand, since M^> and M are solutions of differential equations coming from the 
same system of differential equations (|2.3p . and since they have the same limit as x — > +oo 
for k S: 

Af, M(°) ► J, 



x— >+oo 



it follows that they are related by 

M(x,t,k) = p- 1 (k)D- 1 (x,t)P(k)M^(x,t,k)e^-^ x ' t ^ A ^. (3.12) 

Particularly, at k = m = e~ we have 

/-I 0\ 

A( e f ) = 0, (3.13a) 

\° ° v 

/ l + i -l + i\ 

P" 1 (A ; )D- 1 (x,t)P(fe)|^=M 2(^-1) 2g 2(g-|) . (3.13b) 

Now observe that 

(1 1 l)p- 1 (k)D- 1 (x,t)P(k)\ _ ^ =q(x,t)(l 1 1). 

Combined with ()3.12|) . this implies that the row vector solution /j,(y,t,k) evaluated at k = 
k\ = e"e" takes the value 

/ e -JT(9(^)-l)d£ 

v(y,t,eT) = q ( x ,t)(l 11) 1 

V e JT(<2(^)-i)d^ 

= ( g (x,t)e-r(^)-Dd5 g ( M ) g(x,t)e/. 00 ta«.*)-i)««). (3 . 14) 
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Taking into account that in terms of functions of the variables (x, t) we have u(y, t) = ^jf (y, t) 
(this follows from (12. ip and from the definition (12. lip of the new variable y), the relation (|3.14p 
provides a parametric representation of u. 

Theorem 3.1. Let /x = fj,(y,t,k) = (/ii /i2 ^3) be the solution of the Riemann-Hilbert 
problem (13,7aj) - (13.7cp . where S{k) is the scattering matrix and {v n } are the residues associated 
with the initial data uq(x). 

Then the solution u(x,t) of the Cauchy problem (jl.ip - (jl.2|) for the Degasperis-Procesi 
equation can be expressed in terms of /j,(y,t,k), evaluated at k = e"e" , in parametric form: 

u(y,t) = -lo S ^(y,t,e%), (3.15) 

x(y,t) =y + log^^-(y,t,ef), j = 1 or 2. 

Notice that the structure of the parametric representation (J3.15P is similar to that in the 
case of the Camassa-Holm equation. Moreover, this structure appears also in formulae for 
pure multisoliton solutions given in |21j . 



4. Long-time asymptotics 

The analysis of the long-time behavior of the solution of the IVP is based on the analysis 
of the large-t behavior of the solution of the associated RH problem. The latter can be 
done in the framework of the nonlinear steepest descent method, whose key ingredient is 
the deformation of the original RH problem in accordance with the "signature table" for the 
phase functions involved in the jump matrix. 

The analysis presented in the previous section shows that the structure of the jump matrix 
S, which is 3 x 3, is essentially 2x2: for each straight line of the contour S, a non-trivial 
block of S(k) is 2 x 2, see (|3.4|) . f|3.5[) (under an appropriate change of basis), with only 
one exponential involved for each part, of the form (|3.8|) . Now observe that this exponential 
is essentially the same as in the case of the Camassa-Holm equation (only the constant 2 
is replaced by 3), provided v plays the role of the spectral parameter. Consequently, the 
deformations of the each part of £ are performed in the same way as the deformation of 
the real line in the case of the Camassa-Holm equation, leading to the similar asymptotic 
behavior of the solution, see [Bill]: 

(i) in the sector f > 3, the long-time asymptotics is dominated by the solitons; in the case 
when there are no solitons, the asymptotics in this sector is fast decaying. 

(ii) for < j < 3, the asymptotics has the form of slowly decaying (as t~ 1 ' 2 ) modulated 
oscillations. 

(iii) for — I < j < 0, the asymptotics has the form of a sum of two slowly decaying modulated 
oscillations. 

(iv) j < — |, the asymptotics is fast decaying. 

Besides, there are transition zones between the sectors, where the (slowly decaying) asymp- 
totics is given in terms of the dedicated solutions of the Painleve II equation. 
The details of the asymptotics will be given elsewhere. 
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Figure 3. Sectors in the (x, £)-half-plane with different asymptotics 



References 

[1] R. Beals and R. R. Coifman. Scattering and inverse scattering for first order systems. Coram. Pure Appl. 

Math., 37(l):39-90, 1984. 
[2] Richard Beals, Percy Deift, and Carlos Tomei. Direct and inverse scattering on the line, volume 28 of 

Mathematical Surveys and Monographs. American Mathematical Society, Providence, RI, 1988. 
[3] Anne Boutet de Monvel, Alexander Its, and Dmitry Shepelsky. Painleve-type asymptotics for the 

Camassa-Holm equation. SIAM J. Math. Anal, 42(4):1854-1873, 2010. 
[4] Anne Boutet de Monvel, Aleksey Kostenko, Dmitry Shepelsky, and Gerald Teschl. Long-time asymptotics 

for the Camassa-Holm equation. SIAM J. Math. Anal., 41(4):1559-1588, 2009. 
[5] Anne Boutet de Monvel and Dmitry Shepelsky. Riemann-Hilbert approach for the Camassa-Holm equa- 
tion on the line. C. R. Math. Acad. Set. Pans, 343(10):627-632, 2006. 
[6] Anne Boutet de Monvel and Dmitry Shepelsky. Long-time asymptotics of the Camassa-Holm equation 

on the line. In Integrable systems and random matrices, volume 458 of Contemp. Math., pages 99-116. 

Amer. Math. Soc, Providence, RI, 2008. 
[7] Anne Boutet de Monvel and Dmitry Shepelsky. Riemann-Hilbert problem in the inverse scattering for 

the Camassa-Holm equation on the line. In Probability, geometry and integrable systems, volume 55 of 

Math. Sci. Res. Inst. Publ., pages 53-75. Cambridge Univ. Press, Cambridge, 2008. 
[8] Anne Boutet de Monvel and Dmitry Shepelsky. Long time asymptotics of the Camassa-Holm equation in 

the half-line. Ann. Inst. Fourier (Grenoble), 59(7):3015-3056, 2009. 
[9] Roberto Camassa and Darryl D. Holm. An integrable shallow water equation with peaked solitons. Phys. 

Rev. Lett, 71(11):1661-1664, 1993. 
[10] Roberto Camassa, Darryl D. Holm, and James M. Hyman. A new integrable shallow water equation. 

Hutchinson, John W. (ed.) et al., Advances in Applied Mechanics. Vol. 31. Boston, MA: Academic Press. 

1-33, 1994. 
[11] P. J. Caudrey. The inverse problem for a general n x n spectral equation. Phys. D, 6(l):51-66, 1982/83. 
[12] Adrian Constantin. On the scattering problem for the Camassa-Holm equation. R. Soc. Lond. Proc. Ser. 

A Math. Phys. Eng. Sci., 457(2008) :953-970, 2001. 
[13] Adrian Constantin, Vladimir S. Gerdjikov, and Rossen I. Ivanov. Inverse scattering transform for the 

Camassa-Holm equation. Inverse Problems, 22(6):2197-2207, 2006. 
[14] Adrian Constantin, Rossen I. Ivanov, and Jonatan Lenells. Inverse scattering transform for the Degasperis- 

Procesi equation. Nonlinearity, 23(10):2559-2575, 2010. 
[15] Adrian Constantin and Jonatan Lenells. On the inverse scattering approach to the Camassa-Holm equa- 
tion. J. Nonlinear Math. Phys., 10(3):252-255, 2003. 



16 A. BOUTET DE MONVEL AND D. SHEPELSKY 

[16] A. Degasperis, D. D. Holm, and A. N. I. Hone. A new integrable equation with peakon solutions. Teoret. 

Mat. Fiz., 133(2) T70-183, 2002. 
[17] A. Degasperis and M. Procesi. Asymptotic integrability. In Symmetry and perturbation theory (Rome, 

1998), pages 23-37. World Sci. Publ., River Edge, NJ, 1999. 
[18] P. Deift and X. Zhou. A steepest descent method for oscillatory Riemann-Hilbert problems. Asymptotics 

for the MKdV equation. Ann. of Math. (2), 137(2) :295-368, 1993. 
[19] R. S. Johnson. On solutions of the Camassa-Holm equation. R. Soc. Lond. Proc. Ser. A Math. Phys. Eng. 

Set., 459(2035) T687-1708, 2003. 
[20] Jonatan Lenells. The scattering approach for the Camassa-Holm equation. J. Nonlinear Math. Phys., 

9(4):389-393, 2002. 
[21] Yoshimasa Matsuno. The A^-soliton solution of the Degasperis-Procesi equation. Inverse Problems, 

21(6):2085-2101, 2005. 
[22] Yoshimasa Matsuno. Cusp and loop soliton solutions of short-wave models for the Camassa-Holm and 

Degasperis-Procesi equations. Phys. Lett. A, 359(5) :451-457, 2006. 

*Institut de Mathematiques de Jussieu, Universite Paris Diderot, 175 rue du Chevaleret, 
75013 Paris, France 

E-mail address: aboutet@math.jussieu.fr 

t Mathematical Division, Institute for Low Temperature Physics, 47 Lenin Avenue, 61103 
Kharkiv, Ukraine 

E-mail address: shepelsky@yahoo.com 



